Fermions in the pseudoparticle approach 
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The pseudoparticle approach is a numerical technique to compute path integrals without dis- 
cretizing spacetime. The basic idea is to integrate over those field configurations, which can be 
represented by a sum of a fixed number of localized building blocks (pseudoparticles) . In a couple of 
previous papers we have successfully applied the pseudoparticle approach to pure SU(2) Yang-Mills 
theory. In this work we discuss how to incorporate fermionic fields in the pseudoparticle approach. 
To test our method, we compute the phase diagram of the 1+1-dimensional Gross-Neveu model in 
the large- N limit. 

PACS numbers: 11.15.Tk. 



I. INTRODUCTION 



Recently, there have been a couple of papers propos- 
ing continuum models for SU(2) Yang-Mills theory with 
a small number of physically relevant degrees of free- 
dom. There are ensembles of merons and regular gauge 
instantons [l], there is the pseudoparticle approach 
[H 0, H| and there is a model with calorons with non 
trivial holonomy 0, @] • The basic idea is to restrict the 
Yang-Mills path integral to those gauge field configura- 
tions, which can be represented by a linear superposition 
of a small number of localized building blocks, e.g. instan- 
tons, merons, akyrons or calorons. These models have 
been quite successful when dealing with problems related 
to confinement: the potential between two static charges 
is linear for large separations, there is a confinement de- 
confinement phase transition, and various quantities, e.g. 
the topological susceptibility or the critical temperature, 
are of the right order of magnitude compared to lattice 
results. 

Until now, these models have been applied to SU(2) 
Yang-Mills theory only. In this paper we discuss how to 
include fermionic fields in the pseudoparticle approach. 

The paper is organized as follows. In Sec.[TI]we regular- 
ize the fermionic path integral by considering only those 
field configurations, which can be represented by a lin- 
ear superposition of a fixed number of localized building 
blocks. We point out problems arising in a naive pseu- 
doparticle regularization and propose a solution to these 
problems in form of a slightly different regularization 
scheme. We also discuss possible relations to finite mode 
regularization [1, i, [3. In Sec. DID we test our pseu- 
doparticle method by applying it to a simple interact- 
ing fcrmionic theory, the 1+1-dimensional Gross-Neveu 
model in the large TV-limit [TTJ] . With suitably chosen 
pseudoparticles and after a coupling constant rcnormal- 
ization our pseudoparticle results are in excellent agree- 
ment with analytical results, both for homogeneous chiral 
condensate 12[ 13 1 and for spatially inhomogeneous chi- 
ld, [III]. In Sec. IIV1 we give a summary 



ral condensate 
and a brief outlook regarding the application of the pseu- 
doparticle approach to QCD 



II. FERMIONIC FIELDS IN THE 
PSEUDOPARTICLE APPROACH 

A. Basic principle 

In previous papers @, 0, @ we have discussed in de- 
tail, how to apply the pseudoparticle approach to bosonic 
fields, in particular to the SU(2) gauge field. In this sec- 
tion we propose a method to incorporate fermionic fields. 

The starting point is action and partition function of 
any theory with quadratic fermion interaction: 



S[V>,V,4>] = J dx (^QO)V + 



(1) 
(2) 



where Q is the Dirac operator and <f> denotes any type 
and number of bosonic fields, e.g. the chiral condensate in 
the Gross-Neveu model (cf. Sec. lIII Af or the non-Abelian 
gauge field in QCD. 

To keep close to the spirit of the pseudoparticle ap- 
proach, we regularize the fermionic path integral by con- 
sidering only those field configurations, which can be rep- 
resented by a linear superposition a fixed number of lo- 
calized building blocks: 



j-th pseudoparticle 



3 



(3) 



where rjj are Grassmann valued spinors and Gj are func- 
tions, which are localized in space as well as in time, 
i.e. pseudoparticles. The functional integration over the 
fermionic field configurations is defined via 



Dip Dip. 



(4) 



Moreover, we consider a d+ 1-dimensional periodic space- 
time region of temporal extension Lq and spatial exten- 
sion (Li) d . Note that fermionic fields have to fulfill an- 
tiperiodic boundary conditions in time direction (cf. e.g. 
This implies antipcriodicity for the pseudoparticles 
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Gj, i.e. 

Gj(x + L ,x) = -Gj(x Q ,x). (5) 
After integrating out the fermions one obtains 

•^effective [0] = 

= y*d' 1+1 a ;/:(0)-ln(dct((G 3 |g|Gy))) (6) 

Z (x J £)0e~ SWfe ' !ti ™ M , (7) 

where (Gj \Q\Gj>) is a finite matrix. 

Note that using eigenfunctions of the Dirac operator as 
"pseudoparticles" Gj yields the well known finite mode 
regularization [H, [H, [T(| . 

If det(Q) is real and positive, det(Q) = ^/det(QtQ). 
This suggests another pseudoparticle regularization: 

•^effective [0] = 

= yd d + 1 x£(0)-iln(det((G,|Qtg| Gj ,))). ( 8 ) 

As we will point out in the following, this "Q^Q-regular- 
ization" has significant advantages over the naive "Q-reg- 
ularization" 

B. The Q-regularization versus the 
Q-regularization 

To keep the following arguments as simple as possible, 
we consider all pseudoparticles Gj to be orthonormal, i.e. 
(Gj\Gji) — 5jji. Note that assuming orthonormality is 
not a restriction. Given any set of linearly independent 
pseudoparticles Fj one can easily define suitable linear 
combinations Gj = Mj^Fk, which are orthonormal. Up 
to an additive constant, which is irrelevant for the par- 
tition function, the result for the pseudoparticle regular- 
ized effective action is the same both for pseudoparticles 
Fj and for pseudoparticles Gj: 

ln(dct((G,|Q|G y ))) = 

= ln(det((Aft) fci <F fc |Q|i^)M i / fc /)) = 

= In (dot (VjlQI^-'))) + In (dct . (9) 

—constant 

1. The Problem of the Q -regularization 

In the following we argue that the Q-regularization ([6]) 
is not suited to produce physically meaningful results. 

The problem of the Q-regularization is the following: 
applying the Dirac operator Q to one of the pseudoparti- 
cles Gji in general yields a function, which is (partially) 
outside the pseudoparticle function space span{G n }, i.e. 

QG r {x) = ^, (/i .G,(.r: • hrllyl.n (10) 

k 



with Hj> normalized and Hji _L span{G„}. If 
| J2k a j'kGk\ 3> \hjr\, the situation is uncritical. How- 
ever, as soon as \hj>\ is of the same order of magnitude 
or even larger than \^2 k dj'kGk\ problems arise: when 
computing the matrix elements (Gj\Q\Gji) in ([6]), a sig- 
nificant part of QGj>, hj/Hjt, is simply ignored, just be- 
cause Hj> is perpendicular to the pseudoparticle function 
space span{G„}. 

If Q is hcrmitian, the problem can be made even more 
transparent by diagonalizing the matrix (Gj\Q\Gj'): 

(U f )jk {Gk\Q\Gy) Uk'j' = (Gj\Q\Gj>) = 

= (diag(/ii,^2, ■ ■ ■)) - ( n ) 

v ' or 

where U is a unitary matrix and fXj are the eigenvalues 
of {GjlQlGj,}. Because of dct(WU) = 1, 

det {(GjWGj,)) = HfXj. (12) 

3 

Now the problem is obvious: according to 

H = (Gj\Q\Gj) (13) 

Gj contributes to the pseudoparticle regularized determi- 
nant of Q as eigenmode with eigenvalue fj,j. However, ap- 
plying Q to Gj may yield a function, which is (partially) 
outside the pseudoparticle function space span{G„}, i.e. 

QGj(x) = ^ j Gj(x) + h j H j (x). (14) 

If \hj \ > |/Ltj|, Gj is far from being an eigcnfunction of Q, 
and is, of course, not related to any of the eigenvalues 
of Q. 

The most extreme case is y,j = and hj =/= 0. Gj is then 
an "unphysical pseudoparticle zero mode" : although Gj 
is not a zero mode of Q, and although Q might not even 
have a zero mode, Gj contributes to the pseudoparticle 
regularized determinant of Q with eigenvalue 
fij = 0. Consequently, the effective action © blows up 
to infinity. On the other hand, if we would use additional 
pseudoparticles with non vanishing overlap to QGj, Gj 
would not enter the determinant as zero mode anymore. 

Even if there are no unphysical pseudoparticle zero 
modes, there might still be unphysical low lying pseu- 
doparticle modes, i.e. modes with \hj\ > which also 
spoil numerical results. 

The following example shows that such unphysical low 
lying pseudoparticle modes are quite common. The ex- 
ample in Sec. Ill CI and the ©/-regularized pseudoparti- 
cle Gross-Neveu results from Sec. MI Cl demonstrate that 
these modes usually give rise to wrong and useless results. 

A simple example 

We consider the antihermitian operator d x , 
< x < L G {7 , 8}, antiperiodic boundary conditions. 
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The eigenfunctions of d x are plane waves and the cor- 
responding eigenvalues are given by 



Xr 



2nir 



3 113 

2 ,_ 2' + 2' + 2' 



(15) 



Alternatively, we compute the "pseudoparticle eigen- 
values" fij and \hj\ defined in (fT^)) . where we use over- 
lapping "hat functions" as pseudoparticles. To be more 
precise, we apply antipcriodic B-splinc basis functions of 
degree 2, i.e. 



G 3 {x) 



D 



(2) 

^antipcriodic 



(16) 



j = 0, . . . , L - 1 (cf. Appendix EJ. 

In Table[T] we compare the "true eigenvalues" A r with 
the pseudoparticle eigenvalues We also show \hj\, the 
norm of the overlap of d x Gj to the function space per- 
pendicular to span{G„}. As we have discussed, the pseu- 
doparticle eigenvalues are quite similar to certain true 
eigenvalues as long as ^> \hj\. However, as soon 
as they are of the same order of magnitude, unphysical 
low lying pseudoparticle modes or, in the case of odd L, 
unphysical pseudoparticle zero modes appear. 





I 


, = 7 






I 


. = 8 




r 


A r 


fJ>j 


\hA 


r 






\hA 


1/2 


0.449 i 


±0.449 i 


0.004 


1/2 


0.393 i 


±0.393 i 


0.002 


3/2 


1.346 i 


±1.344 i 


0.135 


3/2 


1.178 i 


±1.177 i 


0.083 


5/2 


2.244 i 


±2.065 i 


1.063 


5/2 


1.963 i 


±1.909 i 


0.599 


7/2 


3.142 i 


0.000 


3.162 


7/2 


2.749 i 


±1.610 i 


2.441 



TABLE I: True eigenvalues X r , pseudoparticle eigenvalues [ij 
and \hi\. 



Note that such unphysical pseudoparticle zero modes 
are not specific for B-splinc pseudoparticles. For odd L 
and any choice of localized, real- valued and uniformly 
distributed pseudoparticles Gj, i.e. 



Gj(x) 



F(x-j)-F(x-j + L), 



(17) 



j = 0, . . . , L - 1, with F* = F and F = for x < and 
x > L, one can easily show that 



G 



zero mo 



do 



(18) 



zero mo 



do IS 



is perpendicular to span{G„}. Therefore, G z 
an unphysical pseudoparticle zero mode of d x 



2. The advantage of the 



The matrix elements {Gj\Q^Q\Gj>) in the G^Q-regu- 
larized effective action © do not suffer from the problem 
discussed in the previous section. The reason is the fol- 
lowing: both the left hand sides (Gj\Q* and the right 
hand sides Q\Gjt) might be outside the pseudoparticle 



function space span{G„}, but they form the same func- 
tion space, span{QG„}, in which their overlap is com- 
puted. 

For example, it is easy to show that any pseudoparti- 
cle zero mode of Q^Q is necessarily a "true zero" mode 
of Q*Q: if G zero mode is a pseudoparticle zero mode of 
Q^Q, then (G zero m ode|Q t Q|G zero mo de) = 0; this implies 

Q|G zcr o mode) = and Q^Q\G zcro mode) = 0. 



3. Another way to motivate the Q -regularization 

It is easy to show that using complete but possibly dif- 
ferent orthonormal sets of functions {F n } and {G„} on 
the left hand side and on the right hand side of any opera- 
tor Q to calculate its determinant yields the determinant 
with exception of a phase factor, i.e. 



det(Q) 



det((Fj|Q|Gj 



(19) 



As we have already stressed, the problem of the Q- 
regularization ^ is that QGy is (partially) perpendic- 
ular to span{G„}. On the other hand, when consider- 
ing operators, where det(Q) is real and positive, we are 
only interested in the absolute value of the determinant. 
Therefore, according to (TT9"|) we propose to consider ma- 
trix elements with different "pseudoparticles" on the left 
hand side and on the right hand side: 

• Right hand side: any choice of pseudoparticles Gy . 

• Left hand side: "pseudoparticles" Fj with 
span{F„} = span{QG„}. 

This assures that Q applied to any pseudoparticle Gy 
yields an clement of the left hand side function space. In 
other words, the best choice for the left hand side is to use 
basis functions Fj with spanji 7 ^} = span{QG n }, because 
other perpendicular/partially perpendicular basis func- 
tions have no/not enough overlap to QGj', and this in 
turn leads to unphysical pseudoparticle zero modes/low 
lying modes. 

One can show that 



det ((F 3 \Q\G 



j' , 



^dct((G,|Qt Q | Gj ,)) ; (20) 



i.e. this line of reasoning also leads to the QtQ-regular- 
ization (T8l) . 



C. The Q^Q- regularization and its relation to finite 
mode regularization 

In the following we point out that there are certain re- 
lations between the G^Q-regularization of the pseudopar- 
ticle approach and finite mode regularization 0, [T^] . 

Let tp n be orthonormalizcd eigenfunctions and X n the 
corresponding eigenvalues of Q^Q, i.e. 



(21) 
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Of course. A„ is real and A„ > 0, because of 



and 



A„ = {ip n \Q ] Q\ip n ) = Q\ipn) 



0. 



(22) H = ( G j\Q Q\ G . 



In the following we consider the eigenvalues ordered ac- 
cording to their absolute value, i.e. |Ao| < |Ai| < ... 
Moreover, we differentiate between "low lying eigenval- 
ues" Ao, • • • , Aa/_i and "large eigenvalues" Am, ■ • - 

• Large eigenvalues 

In applications of the pseudoparticlc approach, 
bosonic fields <p have typical maximum values (cf. 
e.g. [![)■ The reason is both the pseudoparticle reg- 
ularization and the exponential damping factor e~ s 
in the partition function. The same is true for the 
"potential" V(<p) in the Dirac operator 
Q = 7^9^ + V((f)). Therefore, plane waves rje ikx 
with large wave numbers k are approximate eigen- 
functions of Q'Q ("the derivative dominates the 
potential"). Both the cigenfunctions ip n ss r\e %kx 
and the corresponding eigenvalues A„ « ik, 
n = M, . . ., are nearly independent of the bosonic 
fields (p. 

• Low lying eigenvalues 

Low lying eigenvalues A„ and the corresponding 
eigenfunctions ip n , n = 0, . . . , M — 1, may exhibit a 
strong 0-dependence. Because the bosonic fields (p 
have typical maximum values, there is only a finite 
number of low lying eigenvalues. 

Finite mode regularization of det(Q^Q) amounts to 
considering only a finite number of low lying eigenval- 
ues, i.e. 



dct(QTQ) 



N-l 

= " HXj. (23) 

Of course, this equation is not an equality in the usual 
sense. The hope is rather that physical observables can 
be computed correctly via the partition function, when 
using (|23l) in the effective action, and when performing a 
suitable renormalization, 

To exhibit the relation between the C^Q-regularization 
of the pseudoparticle approach and finite mode regular- 
ization, it is convenient to diagonalize the matrix 
{Gj\QiQ\Gjt): 

(tf) jk (G k \QiQ\G k ,)U k/f = (Gj\Q^Q\G f ) = 



= diag(/ii , /i 2 , . . 



33' 



(24) 



where U is a unitary matrix and fij are the eigenvalues 
of {Gj\QtQ\Gj.). Because of Aet{WU) = 1, 



det ((Gj\Q*<Q\Gj>] 



n 



(25) 



In the following we assume that the pseudoparticles Gj 
have been chosen such that the low lying eigenfunctions 
ip n , n = 0, . . . , M — 1, can be approximated. Then 



G, 



^3 



(26) 



A 



3 • 



(27) 



j = 0, . . . , M — 1. That is, the low lying, ^-dependent 
pseudoparticle eigenfunctions Gj and the pseudoparticle 
eigenvalues fj,j are nearly identical to the true eigenfunc- 
tions ijjj and eigenvalues Xj. Therefore, the contribution 
of these pseudoparticle eigenmodes to the pseudoparticle 
regularized determinant of Q^Q is identical to the con- 
tribution of the corresponding true eigenmodes in finite 
mode regularization (cf. (|2"3"|) and ([23)) ). 

All other pseudoparticlc eigenfunctions are essentially 
linear combinations of the remaining true eigenfunctions: 



G, 



n=M 



n=M 



1, 



(28) 



j = M, . . . (there is no contribution from ipo,..., ipM-i> 
because the pseudoparticle eigenmodes are orthonormal 
due to the hermiticity of {Gj\Q^Q\Gj>)). The corre- 
sponding pseudoparticle eigenvalues fij are given by 



= (Gj\Q^Q\Gj) 



\ C jn 



A„ > Am, (29) 



j = M,..., where and 
are two important points: 



have been used. There 



• Gj and fij, j — . . ., arc approximately indepen- 
dent of <p. That is, although the difference between 
these pseudoparticlc eigenvalues fij and the true 
eigenvalues Xj, j = M,. . ., might be large, the pseu- 
doparticle effective action merely differs by an ad- 
ditive constant, when compared to finite mode reg- 
ularization. Such an additive constant is, of course, 
irrelevant for the partition function. 

• fJ-j > Am, j = M,..., i.e. fij is large. Therefore, 
possibly present weak ^-dependencies of Cj n or Xj 
or weak contributions from low lying eigenmodes 
ip n , n = 0, . . . , M — 1, do not have a strong impact 
on fij. Assuming that there are not too many of 
these large pseudoparticlc eigenmodes, their (^-de- 
pendence is essentially negligible in the pseudopar- 
ticle regularized determinant of Q^Q. 

It is also instructive to discuss the shortcoming of the 
Q-regularization in this context. For hermitian Q the 
above arguments and equations are quite similar: just 
replace Q^Q by Q in to (gSD with exception of (|2"2"j) . 
which is of course not valid. (|29|) must be replaced by 



{Gj\Q\Gj 



E 

n=M 



A,, 



(30) 



The big difference compared to the Q^Q- regularization 
is that the true eigenvalues Xj are, in general, not posi- 
tive. Therefore, the corresponding pseudoparticle eigen- 
values fij, j — M, . . ., are not necessarily large: due to 
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intricate cancellations between positive and negative A,,-, 
there may very well be small pseudoparticle eigenvalues 
fij, which are strongly affected by small changes of Cj n 
or Xj , because of a change of <j>. 



Another simple example 

In this example we compare the Q^Q-regularization, 
the Q-regularization and finite mode regularization. To 
this end, we compute 



ln(det(d x + o)) - L = 7 



ln^det(0 x + trJJ, (31) 
a real, < x < L — 7, antiperiodic boundary conditions. 

• Finite mode regularization 

The eigcnfunctions of d x + a are plane waves and 
the corresponding eigenvalues are given by 



2nir 



(32) 



According to ([2"3"]l the result is 
In (det (d x = 



N-l 



/ 27r(n+l/2) 



(33) 



when using the 2N "lowest lying eigenmodes" . For 
N = 3 and N = 4 it is shown as a function of a in 
Fig.ffl 

• Q^Q-regularization and Q-regularization 

We use orthonormalized B-splinc pscudoparticles of 
degree 2 as pseudoparticles Gj. The Q^Q-result, 

\ In ( det ({Gj | (d x + a) ' (d x + a) \G r ))), (34) 
and the Q-result, 

ln(dct((G J |(9 ;E + a)|G / ))), (35) 

are shown as functions of a in Fig. [TJ 

The results obtained with Q^Q- regularization and with 
finite mode regularization arc nearly identical with excep- 
tion of a shift along the vertical axis. This shift is due 
to the different number of degrees of freedom: 7 for the 
pseudoparticle regularization, 6 and 8 for finite mode reg- 
ularization. The Q-regularized result, on the other hand, 
has a completely different shape, especially for small val- 
ues of a. As explained above, the reason are unphysical 
low lying pseudoparticle modes. At a = an unphysical 
pseudoparticle zero mode causes a singularity. 



D 2 - 




finite mode regularization, N = 3 and N = 4 
Q f Q-regularization 
Q-regularization 



0.2 



0.4 



0.6 



0.! 



FIG. 1: L — 7. k^det^ + a)) as function of a computed 
via finite mode regularization, Q-regularization and Q-reg- 
ularization. 



Computing approximate eigenfunctions via the 
pseudoparticle approach 

The pseudoparticle approach can also be used to com- 
pute approximate eigenfunctions of Q*Q (if Q is hermi- 
tian, one can also use Q instead of Q^Q): 

• Let Gj be a set of orthonormalized pseudoparticles. 

• Compute the matrix {Gj\Q^Q\Gji) and diagonalize 
it according to (f2"4"| . 

• Compute \hj\, defined by 

Q^QGjix) = HjGj(x) + hjHj(x) (36) 
with Hj normalized and Hj _L span{G n }. 

• If \fij \ 3> \hj\, the pseudoparticle cigenfunction Gj 
is close to a true eigenf unction of Q^Q and the cor- 
responding pseudoparticle eigenvalue fj,j is close to 
a true eigenvalue. 

• If \fij\ ~ \hj\ or if \fXj\ < \hj\, the pseudoparticle 
eigenfunction Gj is not close to a true eigenfunction 
of QtQ. 



The simple example continued 



We continue the example of Sec. Ill B 11 by computing 
the pseudoparticle eigenfunctions Gj of the antihermitian 
operator d x for L = 7. 

In Fig. [2] we compare both the real parts and the imag- 
inary parts of the pseudoparticle eigenfunctions and the 
true eigenfunctions (l/VL)e XrX with A r given by (fT^j) . 
As long as |^-| > \hj\ (r = 1/2 and r = 3/2; cf. Table[Q) 
pseudoparticle results and analytical results are essen- 
tially identical. However, if \/j,j\ < \hj\ (r — 7/2), Gj is 
not close to a true eigenfunction of d x anymore. 
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r = 3/2 



N 




r = 5/2 



= 7/2 





FIG. 2: L = 7. Real parts and imaginary parts of the pseu- 
doparticle eigenfunctions and the true eigenfunctions of d x as 
functions of x. 



III. THE 1+1-DIMENSIONAL GROSS-NEVEU 
MODEL IN THE PSEUDOPARTICLE 
APPROACH 

In the following, we test the pseudoparticlc approach 
by considering a simple interacting fermionic theory, 
the 1+1-dimensional Gross-Neveu model in the large-TV 
limit. We compute the phase diagram both for homoge- 
neous and for spatially inhomogencous chiral condensate. 



A. The l + l-dimensional Gross-Neveu model in 
Euclidean spacetime 

The Gross-Neveu model is a four fermion interact- 
ing theory with N degenerate flavors. Action and parti- 
tion function are given by 



S = 



r ( N 

J d 2 x ( ^ (ri) ( 7 o(c>o + M)+7i3i)?/> ( 



N 



N 




(37) 
(38) 



where N is the number of flavors, [i is the chemical poten- 
tial and g is the coupling constant. The Dirac matrices 
fulfill {7^,7^} = 25^, e.g. 7 = cr 1 and 7! = a 3 . 

To get rid of the four fermion interaction, one usually 
introduces a real scalar field a: 



S' 



n=l 



-Q (Dirac operator) 



Z oc 



Integrating out the fermions yields 

^effective = N j (fx O 2 - hi ( det(Q) 



Z oc 



(39) 
(40) 

(41) 
(42) 



where A = Ng 2 . 

From now on we consider the large- N limit, i.e. 
N — > oo and A = constant. Note that due to 
^effective N, only a single field configuration con- 
tributes to the partition function (|4"2"|) . It can be de- 
termined by minimizing -^effective with respect to a. 

Moreover, one can show that in the large- N limit 



iV 



(n) 



(43) 



i.e. the scalar field a is proportional to the chiral conden- 
sate. 



B. B-spline pseudoparticles 

For the following computations we use a large number 
of overlapping "hat functions" as pseudoparticles: we ap- 
ply products of antipcriodic and periodic B-spline basis 
functions of degree 2, i.e. 

Gj Ji(«0,a?l) = ^ojantiperiodic^o)^^^^^!) (44) 

j = 0, . . . ,L - 1, ji = 0, . . . , L x - 1 (cf. Appendix EJ), 
where Lq x Li is the extension of the periodic spacetime 
region (Lq and L\ arc chosen to be integers). Fig.[3]shows 
the "B-spline pseudoparticlc" Goo- 




FIG. 3: B-spline pseudoparticle Go.o as a function of (a;o,a;i). 
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Why using B-spline pseudoparticles? 

• The intention of this section is to test, whether the 
pseudoparticle approach as presented in the previ- 
ous section is suited to deal with fermionic fields. 
Therefore, we need pseudoparticles, which form a 
"sensible set of field configurations" , i.e. pseudopar- 
ticles, which are able to approximate any not too 
heavily oscillating field configuration. For B-splinc 
pseudoparticles, i.e. for a piecewise polynomial ba- 
sis, this is certainly the case. Note that it is not 
our intention to determine a small number of phys- 
ically important fermionic field configurations; this 
will be part of an upcoming paper, where we apply 
the pseudoparticle approach to QCD. 

• Since B-spline pseudoparticles are piecewise poly- 
nomial functions, the matrix elements 
(Gj ,ji\Q^Q\Gj' , j{) can be calculated analytically. 

• B-spline pseudoparticles of degree 2 with uniform 
knot vectors tj = j are localized within a space- 
time region of extension 3x3. Therefore, when 
using a large number of pseudoparticles the ma- 
trix (GjgjilQ'QlGj' j') is sparse, which is benefi- 
cial from a numerical point of view. 



C. Homogeneous chiral condensate 

In this section we consider a homogeneous chiral con- 
densate, i.e. a = constant. 

The Q^Q-regularized pseudoparticle effective action is 
given by 



0.0, L n x L, = 8 x 1 44 
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effective 
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-ln(det((G JOJ1 |QtQ| G 



3q >3i 



(45) 



To determine the chiral condensate a for given temper- 
ature T = I/Lq and chemical potential /x, one has to 
minimize this expression with respect to a. 

Of course, numerical results strongly depend on the 
number of pseudoparticles applied. To extract physically 
meaningful results, a coupling constant renormalization 
is necessary, i.e. we have to choose A in accordance with 
the number of degrees of freedom. A possible way of do- 
ing that is to perform finite temperature computations 
with temporal extension Lq = 8 much smaller than spa- 
tial extension L\ = 144 at \i = 0.0 for different values of 
A. The resulting a as a function of A is shown in Fig. [4ji. 
For all further computations we fix the coupling constant 
A at that value, where the chiral condensate a just van- 
ishes: A = 1.153. The scale is now set, i.e. 1/Lq = 1/8 
corresponds to the critical temperature of chiral symme- 
try breaking at fi = 0.0. 

We can now compute the chiral condensate a at "arbi- 
trary" temperature T = l/L$ (Lq integer) and chemical 
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FIG. 4: Q+Q-regularization, Li = 144. a) fj, = 0.0, L = 8. 
a as a function of A. b) A = 1.153. The pseudoparticle phase 
diagram for homogeneous chiral condensate together with the 
exact phase boundary. 



potential fi, to determine the phase diagram of the Gross- 
Neveu model. As usual, we express T and \x in terms of 
do, which is the T = 0.0 and n = 0.0 value of the chi- 
ral condensate a (to determine <7o we have performed a 
"zero temperature computation" with large temporal ex- 
tension Lq = 48). The resulting pseudoparticle phase 
dia gram is in excellent agreement with analytical results 



T2M (cf- Fig.Hb). 



Another possibility to adjust the temperature T/cto is 
to change the coupling constant A, while the extension 
of the spacetime region Lq x Li is kept constant. This 
is similar to what is usually done in lattice calculations. 
The scale is then set via <7o, which is A-dcpcndcnt. Using 
this method the resulting pseudoparticle phase diagram 
is also in excellent agreement with analytical results. 

In Sec. Ill Bl we have pointed out that the Q-regulariza- 
tion is not suited to produce physically meaningful re- 
sults. To demonstrate that this is indeed the case, we 
perform similar computations with the Q-regularized ver- 
sion of the effective action, 
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(46) 



For odd Lo unphysical pscudoparticle zero modes ren- 
der the results completely useless: there is no chirally 
symmetric phase and, therefore, no sensible phase dia- 
gram. For even Lq there are no unphysical pseudoparticle 
zero modes, but still unphysical low lying pscudoparticle 
modes. The resulting phase diagram is not in quantita- 
tive agreement with analytical results [H,[l3| (cf. Fig. [5]). 
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0.6 
0.5 
0.4 
0.3 
0.2 
0.1 

o 



a ^ 



exact phase boundary 
pseudoparticle results 



G = 



0.2 0.4 0.6 0.8 1 1.2 1.4 
chemical potential u/a 



1.6 



FIG. 5: Q-regularization, A = 0.894, L even, L\ = 144. The 
pseudoparticle phase diagram for homogeneous chiral conden- 
sate together with the exact phase boundary. 



D. Spatially inhomogeneous chiral condensate 

In this section wc consider a spatially inhomogeneous 
chiral condensate, i.e. a = <r(xi). This gives rise to a new 
so called crystal phase. In this phase the chiral conden- 
sate is given by 



cr(xi) 



A, 



2 sn(Ax, k)cii(Ax, k) 
dn(^4x, k) 



(47) 



where sn, cn and dn are Jacobi elliptic functions and A 
and k are functions of T and fi fl5| . 



1. The right phase boundary of the crystal phase 

In addition to the fermionic fields we also represent the 
chiral condensate in terms of B-splinc pseudoparticles: 



cr(xi) 
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with uniform knot vector tj = 3j and L\ a multiple of 3 
(cf. Appendix [A"|l . Note that we consider three times as 
many degrees of freedom in Xi-direction for the fermionic 
fields than for the chiral condensate. The reason is the 
following: in general, even low lying eigenfunctions of 



Q^Q exhibit more oscillations than the chiral condensate 
cr; on the other hand, the fermionic pseudoparticles must 
be able to approximate these low lying eigenfunctions, 
in order to produce correct results; therefore the repre- 
sentation of the fermionic fields must be finer than the 
representation of the chiral condensate. 

In accordance with Sec. IIII Cl we choose A = 1.153 and 
L\ = 144. We determine the right phase boundary of the 
crystal phase by computing the eigenvalues of the Hessian 
matrix of the Q^Q-regularized effective action (|4"5"j) with 
respect to a n at a = 0, i.e. 



H,, 



_d d_ 



ffcclr 



(49) 



Negative eigenvalues indicate inhomogeneous perturba- 
tions of the chiral condensate a = 0, which decrease the 
effective action. This in turn is a clear sign of a crystal 
phase. Of course, the argument does not work the other 
way round: if there are no negative eigenvalues, the chi- 
ral condensate is not necessarily vanishing; a — may as 
well be a local minimum of the effective action; however, 
a recent investigation of the Gross-Neveu model on the 
lattice strongly suggests that this is not the case (l7j . 

A minor source of error when computing the right 
phase boundary is the finite extension of the space di- 
mension. The problem is that the period of the analyti- 
cally obtained chiral condensate (|47p might not "fit in the 
periodic xi-direction" . If the period of the chiral conden- 
sate is a multiple of the spatial extension L\, everything 
works fine. However, if it is roughly halfway between two 
multiples, results differ a little bit (cf. Fig. where the 
lowest eigenvalue of H is shown as a function of /j/uq 
[Lq = 36, i.e. T/cto = 0.126]; for a space dimension of 
infinite extension the graph would be smooth, instead of 
exhibiting certain periodic oscillations). Note that the 
same behavior has been observed on the lattice [IT} . 

To get rid of these finite size effects, we fit a parabola 
such that it just touches the curve of the lowest eigen- 
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FIG. 6: A = 1.153, L x Li = 36 x 144, i.e. T/a = 0.126. 
The lowest eigenvalue of the Hessian matrix H as a function 
of n/a . 
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FIG. 7: QtQ-regularization, A = 1.153, L x = 144. The pseu- 
doparticle phase diagram for spatially inhomogeneous chiral 
condensate together with the exact phase boundary. 



FIG. 8: A = 1.153, fj,/<r = 0.704, L x L 1 = 24 x 144, i.e. 
T/ao = 0.19. Pseudoparticle chiral condensate and the corre- 
sponding analytically obtained chiral condensate as functions 
of Xl. 



values from below (cf. Fig. [5]). The root of this parabola 
is then taken as the corresponding (^i/oo)- value of the 
right phase boundary. The resulting right phase bound- 
ary is shown in Fig. [JJ together with the analytical result 
[15( 1 . Up to /i/(To ~ 1.5 there is excellent agreement. For 
larger values of fi/a cutoff effects give rise to certain 
deviations: the analytically obtained chiral condensate 
(|4"T)) oscillates heavily and, therefore, cannot be repre- 
sented by Li/3 = 48 B-spline pseudoparticles anymore. 
Of course, applying a larger number of pseudoparticles 
in xi-direction allows to extract correct results for larger 
values of /i/orj. 

2. The left phase boundary of the crystal phase 

Proceeding in the same way to determine the left phase 
boundary docs not yield correct results. One merely ob- 
tains the phase boundary of the phase diagram for ho- 
mogeneous chiral condensate (cf. Fig.[4i)). The reason is 
that minima of the effective action for a = constant ^ 
are also minima of the effective action for varying a (how- 
ever, not necessarily global minima), i.e. the Hessian ma- 
trix at such values for a is positive definite. To get the 
correct phase boundary, one has to perform a minimiza- 
tion of the pseudoparticle effective action (|45|) with re- 
spect to er„. Again this is i n ag reement with what has 
been observed on the lattice [l7[ • 

The resulting left phase boundary is shown in Fig. [TJ 
It is in excellent agreement with the analytical result [l5| • 



shown in Fig. M W^o = 0.704, T/a Q = 0.19). There is 
excellent agreement between pseudoparticle and analyti- 
cal results. 

It is also interesting to compute the chiral conden- 
sate for different values of fx/cro along a line of constant 
T/ao- Results are shown in Fig. [S] as functions of x\ for 
T/a = 0.141. In agreement with [3, EH, tne chiral 
condensate changes from a sin-like behavior inside the 
crystal phase (cf. Fig. [5)3) to a kink-antikink structure, 
when approaching the left phase boundary (cf. Fig.[5K)- 
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FIG. 9: A = 1.153, L x Li = 32 x 144, i.e. T/a = 0.141. 
a as a function of xi. a) Close to the left phase boundary: 
fj>/ao = 0.636. b) Inside the crystal phase: fi/cro = 0.660. 



3. The chiral condensate 



Another check of the pseudoparticle approach is to 
compare the chiral condensate obtained by minimizing 
the Q^Q-regularized effective action ([4*5)l with the ana- 
lytically obtained chiral condensate (JUJ- An example is 



IV. CONCLUSIONS AND OUTLOOK 

In this paper we have presented a method to incorpo- 
rate fcrmions in the pseudoparticle approach. 
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We have pointed out that a naive pseudoparticle rcgu- 
larization, the Q-regularization, is not suited for produc- 
ing useful numerical results. The problem of this regular- 
ization is that applying the Dirac operator to the pseu- 
doparticles yields functions, which are partially outside 
the pseudoparticle function space. This gives rise to un- 
physical low lying modes or, in extreme cases, unphysical 
pseudoparticle zero modes, which spoil the determinant 
in the effective action. 

A slightly different regularization scheme, the Q^Q- 
regularization, does not suffer from this problem. To 
test this pseudoparticle regularization, we have computed 
the phase diagram of the 1+1-dimcnsional Gross-Neveu 
model in the large- N limit, both for homogeneous and 
for spatially inhomogeneous chiral condensate. Our pseu- 
doparticle results are in quantitative agreement with ana- 
lytical results. We have given a couple of arguments indi- 
cating certain relations between the Q'Q- regularization 
and finite mode regularization, which might explain these 
excellent results. 

The next step is to apply the pseudoparticle approach 
to QCD. Of course, instead of using a large number of 
overlapping hat functions as pseudoparticles, as it has 
been done in this paper, the goal will rather be to use a 
small number of physically relevant fermionic pseudopar- 
ticles, probably pseudoparticles, which have significant 
overlap to low lying eigenfunctions of the QCD Dirac op- 
erator. 

Current research includes a study of chiral symmetry 
breaking in QCD in the pseudoparticle approach. We do 
this by computing the low lying eigenvalues of the Dirac 
operator at different temperatures and relate the results 
to the chiral condensate via the Banks-Casher relation. 
Of course, the goal is to obtain a model which exhibits 
chiral symmetry breaking and a confinement deconfinc- 
ment phase transition at the same time. Such a model 
could be useful for computing obscrvables, which are dif- 
ficult to access in lattice calculations, e.g. pion masses or 
decay constants with realistically light quark masses. 
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= { I otherwise t '- X< < A1 > 

t J+k - tj 

/ W1 7 fifths) (A2) 

with suitably chosen knot vector ... < t j < tj +1 < . . ., 
e.g. /, j (cf. e.g. \EM l 

For < x < L, L integer, antiperiodic and periodic 
B-spline basis functions with uniform knot vector tj = j 
are given by 

4antipenodi>) = B?\x) - B?} L {x) (A3) 

^5 oriodic (z) = B?\x)+Bf} t {x), (A4) 

j = 0, . . . , L — 1. For L = 6 and k = 2 they are shown in 
Fig.rini 

a) antiperiodic B-spline basis functions, degree 2 - L = 6 
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b) periodic B-spline basis functions, degree 2 - L = 6 
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APPENDIX A: B-SPLINE BASIS FUNCTIONS 

B-spline basis functions are piecewise polynomial func- 
tions of degree fc, which are C fc_ ^continuous: 



FIG. 10: L — 6. a) Antiperiodic B-spline basis functions of 
degree 2. b) Periodic B-spline basis functions of degree 2. 
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